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Diagnostics  informed  by  observations  and  compared  with  model 
output 


Current  meters  (Global  Multi-Archive  Current  Meter 
Database; 

http://stockage.univ-brest.fr/scott/GMACMD/updates.html) 


9  Mean  vertical  structure  of  kinetic  energy 
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9  Mean  vertical  structure  of  kinetic  energy 


Satellite  altimetry  (Archiving,  Validation  and  Interpretation 
of  Satellite  Oceanographic; 

http://www.aviso.oceanobs.com/es/data/index.html) 


9  Surface  kinetic  energy 

9  Eddy  length  scales  (inverse  first  centroid  of  kinetic  energy 
power  spectrum) 

9  Sea  surface  height  variance 
9  Intensified  jet  positions  (via  Kelly  et  at.,  2007) 
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Wave  drag  scheme  choices 

1)  What  is  the  range  of  wavenumbers  over  which  the  internal 
waves  are  not  evanescant? 


f/U  ~  1 0-4m-1  <  \k\  <  N/U  ~  1  (T1  m"1 . 

Here, 

o  f  is  the  Coriolis  parameter 
o  N  is  the  buoyancy  frequency 
o  U  is  the  velocity  near  the  seafloor 
®  \k\  is  the  wavenumber  of  the  internal  wave 
Scott  et  at.  (201 1 )  used  a  range  that  went  down  to  f/U  ~  1 0“ 
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2)  Which  wave  drag  parameterizations  are  there  to  choose  from? 


Using  a  momentum  sink: 


9  Implement  in  wavenumber  space;  e.g.,  Bell  (1975) 
9  Implement  in  physical  space;  e.g.,  Gamer  (2005) 
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Using  a  momentum  sink: 


9  Implement  in  wavenumber  space;  e.g.,  Bell  (1975) 
9  Implement  in  physical  space;  e.g.,  Gamer  (2005) 


Features  of  Gamer  (2005)  vs  those  of  Bell  (1975) 


9  Garner  (2005)  -  allows  for  topographic  blocking,  but  does 
not  depend  on  Coriolis 

9  Bell  (1975)  -  does  not  allow  for  topographic  blocking,  but 
does  depend  on  Coriolis 

9  Both  schemes  -  depend  on  stratification,  velocity,  and 
underlying  topographic  features  and  assume  f  <  N 
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2)  (cont...)  Comparison  of  the  Bell  (1975)  and  Gamer  (2005) 
schemes 


We  choose  to  use  the  Gamer  (2005)  scheme,  but  the  Bell 
(1975)  scheme  yields  similar  results  (offline) 


a)  (Garner  scheme,  N  from  WOA)  b)  (Garner  scheme,  N  from  HYCOM) 


'-6 


Introduction 

oooooooo 


Putting  wave  drag  into  an  ocean  model 
0000*00 


Energy  budget 
oooooo 


Model  evaluation 
oo 


Summary 


Wave  drag  scheme  choices 


3)  Where  do  we  apply  wave  drag? 


»  Is  the  model  numerically  stable  when  wave  drag  is  applied 
everywhere? 

®  Is  it  possible  and  does  it  make  sense  to  apply  wave  drag 
everywhere? 

Interpolate  over  topographic  slopes  that  are  supercritical? 
Apply  wave  drag  only  in  abyssal  hill  regions?  Apply  wave  drag 
only  in  regions  deeper  than  500  meters?  . . . 
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Wave  drag  scheme  choices 

4)  Estimate  the  input  parameters  for  the  wave  drag  scheme  of  your 
choice 


9  Integrate  Goff  and  Jordan  (1988)  abyssal  hill  power  spectrum,  weighted  by 
wavenumbers  from  (1) 

9  parameters  for  power  spectrum  from  Goff  and  Arbic  (201 0)  and  Goff  (201 0)  in 
abyssal  hill  regions 

9  use  a  machine  learning  algorithm  (Wood,  2006)  to  fill  in  the  non-abyssal  hill 
regions 


c)  10x(2,1)  component  of  info,  tensor  d)  (2,2)  component  of  info,  tensorf  -2.5 


180°W  108°W  36°W  36°E  108°E  180°W  108GW  36°W  36°E  108°E  180°W 


a)  (1,1)  component  of  info,  tensor 
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Wave  drag  scheme  choices 

5)  How  should  the  momentum  be  deposited  vertically? 


®  Is  there  observational  evidence  for  enhanced  turbulence,  if 
not  lee  wave  drag,  in  the  bottom,  say,  500  meters?  (see 
Naveira-Garabato  etal.,  2012) 

®  Is  there  evidence  that  there  needs  to  be  a 
depth-dependent  vertical  deposition  of  momentum? 

(Polzin  (2009)  suggests  that  there  is  and  the  Gamer 
(2005)  scheme  is  capable  of  doing  this) 

®  Are  there  locations  where  a  non-trivial  vertical  deposition 
of  momentum  is  important?  (will  not  be  addressed  here) 
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Multiply  the  momentum  equations  by  p  and  take  a  dot  product 
with  velocity,  u;  then  integrate  over  the  globe 
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Multiply  the  momentum  equations  by  p  and  take  a  dot  product 
with  velocity,  u;  then  integrate  over  the  globe 


P EKtime  +  P EKadvection  —  P pressure  +  P input  ~  P output  +  Cek 


(3) 
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Global  Integrals  of  Input/Output  Terms  in  TW=  1012W 


Pek  time  +  P EKadvection  —  Cek^Ep  +  P pressure 
+  P Wind  -  PbD  ~  PwD  ~  Pw  -  PhV 


(4) 


WD? 

Wind 

Buoy 

BD 

WD 

VV 

HV 

no 

0.87 

0.066 

0.31 

N/A 

0.29 

0.29 

yes 

0.87 

0.066 

0.14 

0.40 

0.28 

0.26 

Inputs  vs  Outputs: 

»  5%  imbalance  (outputs  less  than  inputs)  without  wave  drag 
«  15%  imbalance  (inputs  less  than  outputs)  with  wave  drag 
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Mass  conservation  equation  potential  energy  equation 
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Mechanical  energy  budget  from  the  continuity  and  momentum  equations 

Global  Integrals  of  Mechanical  Energy  Budget  Terms  in 

TW  1012W 


PEKtime  +  P Eptime  +  PEKadvection  +  PEPadvection  —  (8) 

P pressure  +  P diffusive  +  P input  ~  P output  +  Ce^Ep 


KEadv. 

PEadv. 

press . 

diffuse 

Ei^EP 

input 

output 

-.00284 

.174 

<  .001 

.00309 

.0865 

.868 

1.06 

7%  imbalance  of  mechanical  energy  budget  we  ignore: 
«  partial  time  derivatives  of  KE  and  PE 
9  along-isopycnal  contributions  to  power  associated  with 
buoyancy  diffusion 
9  compressibility 
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Taylor  diagrams  of  all  five  diagnostics 

Does  wave  drag  ever  make  the  model  simulations  in  worse 
agreement  with  diagnostics  informed  by  observations? 


Observations,  1/12°  HYCOM  without  wave  drag,  1/12°  HYCOM  with  wave  drag,  1  /25°  HYCOM  without  wave 
drag,  1  725°  HYCOM  with  wave  drag,  1  710°  POP  without  wave  drag  ( Taylor ,  2001) 


a)  Kinetic  energy  profiles  Taylor  diagram  b)  Intensified  jet  positions  Taylor  diagram  c)  SSH  variance  Taylor  diagram 
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There  are  several  ways  wave  drag  impacts  the  model 

®  active  feedback  on  velocities  and  stratification 

9  diapycnal  diffusivity  is  generally  enhanced  all  the  way  up  to 
the  surface 

9  substantially  less  bottom  drag  dissipation  with  wave  drag, 
and  wave  drag  cannot  be  substituted  for  by  boosting 
bottom  drag 

9  all  other  mechanical  energy  budget  terms  are  spatially 
altered,  but  changed  by  little  in  their  global  integrals 

9  wave  drag  either  improves  the  model  or  does  not  make  the 
model  worse 


Non-input/output  mechanical  energy  budget  terms 


a)  Pressure  work  [W  m 2]  b)  Buoyancy  diffusion  [W  m 2]  c)  I E  to  PE  conversion  [W  rrr] 


d)  PE  advection  [W  m'2]  e)  KE  advection  [W  m'2J 
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SST  bias 


a|  SST  percent  difference  1112th  HYCOM  without  wave  drag  b  SST  percent  difference  1125th  HYCOM  without  wave  drag  c)  SST  percent  difference  l/ioth  POP  without  wave  drag 
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d]  SST  percent  difference  1/121  HYCOM  with  wave  drag  e  SST  percent  difference  1125th  HYCOM  with  wave  drag 
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